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. ABSTRACT 
O ■ 

^: 

Ph, We present results of an extensive mapping survey of N2H'^(l-0) in about 

O I 60 low mass cloud cores already mapped in the NH3(1,1) inversion transition 

^ I line. The survey has been carried out at the FCRAO antenna with an angular 

I resolution of 54", about 1.5 times finer than the previous ammonia observations 

I made at the Haystack telescope. The comparison between N2H+ and NH3 maps 

^ I shows strong similarities in the size and morphology of the two molecular species 

I indicating that they are tracing the same material, especially in starless cores. 

Cores with stars typically have map sizes about a factor of two smaller for N2H"'' 
than for NH3, indicating the presence of denser and more centrally concentrated 
gas compared to starless cores. The mean aspect ratio is ~ 2. Significant 
correlations are found between NH3 and N2H^ column densities and excitation 
temperatures in starless cores, but not in cores with stars, suggesting a different 
chemical evolution of the two species. Starless cores are less massive (< M^ir > 
~ 3 Mq) than cores with stars (< M^ij. > ~ 9 Mq). Velocity gradients range 
between 0.5 and 6 km/s/pc, similar to what has been found with NH3 data, and 
the ratio (3 of rotational kinetic energy to gravitational energy have magnitudes 
between ~ 10~^ and 0.07, indicating that rotation is not energetically dominant 
in the support of the cores. "Local" velocity gradients show significant variation 
in both magnitude and direction, suggesting the presence of complex motions 



- 2 - 



not interpretable as simple solid body rotation. Integrated intensity profiles of 
starless cores present a "central flattening" and are consistent with a spherically 
symmetric density law n oc r~" where a = 1.2 for r < rtreak and a = 2 for r > 
^break, where Tbreak ~ 0.03 pc. Cores with stars are better modelled with single 
density power laws with a > 2, in agreement with observations of sub millimeter 
continuum emission. Line widths change across the core but we did not find 
a general trend: there are cores with significant positive as well as negative 
linear correlations between Av and the impact parameter b. The deviation in 
line width correlates with the mean line width, suggesting that the line of sight 
contains ~ 10 coherence lengths. The corresponding value of the coherence 
length, ~ 0.01 pc, is similar to the expected cutoff wavelength for MHD waves. 
This similarity may account for the increased "coherence" of line widths on 
small scales. Despite of the finer angular resolution, the majority of N2H"'" 
and NH3 maps show a similar "simple" structure, with single peaks and no 
elongation. 



Subject headings: molecular data - ISM: clouds, molecules, structure - radio 
lines: ISM 



1. Introduction 



Dense cores in dark clouds have been extensively studied through observations of 
the inversion transition lines of ammonia ( [Myers fc Benson 1983| ; [Benson fc Myers 1989 



hereafter BM89) and other high density tracers, including CS ( [Zhou et al. 198^ , C2S 
( [Suzuki et al. 1992[ ) and HC3N ( [b'uUer fc Myers 1993[ ). Molecular emission maps have 
improved our understanding of cloud structure and have given us insights on the initial 
conditions of the star forming process. It is now well established that low mass cores as 
mapped in NH3 lines are about 0.1 pc in size, have kinetic temperatures of about 10 K, 
and gas number density of ~ 3x10^ cm~^ (BM89). Generally, cores have elongated maps 
with typical aspect ratio of 2 ( Myers et al. 199"l[ ; [Ryden 1996 ). Goodman et al. (1993) used 



ammonia maps to measure velocity gradients; they found typical magnitudes between 0.3 
and 4 km s 
of cores. 



^ pc ^ and conclude that rotation is not energetically dominant in the support 



The upgrade of the 37-m Haystack telescope Q at 3 mm (Parvainis fc Salah 1994[ ) 



^ Radio Astronomy observations at the Haystack Observatory of the Northeast Radio Observatory 
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enabled high spatial (and spectral) resolution observations of the molecular ion N2H+ 
towards those cores already studied in NH3 ( [Benson, Caselli fc Myers 199^ , hereafter 
BCM98). The J = 1 ^ rotational transition of diazenylium has been detected in most 
(94%) of the cores indicating that this species is widespread and easy to detect. Moreover, 
a good correlation between N2H"'" and NH3 velocities and line widths indicates that the two 
species are probably tracing the same material. N2H"'" is known to be a selective tracer 
of quiescent gas ([Turner fc Thaddeus 197"7| ; [Womack et al. 1992[ ; [Bachiller 1996[) and is 
particularly suitable for studying the structure and kinematics of cold star forming cores. 
Being an ion, diazenylium is also important to trace the ionized gas and to give information 
about the coupling between ions and neutrals in star-forming dense cores (BCM98). We 
point out that N2H+, formed through the ion-molecule reaction N2 + H3' and mainly 
destroyed by CO and electrons, traces molecular nitrogen, N2, which is a major repository 
of nitrogen. N2H+ is thought to be a "late depleter" ( [Bergin fc Langer 1997[ ; [Aikawa et al 
'20011 [Uaselh et al. 200r5[) and so is a good tracer of dense core gas. 



In this paper we present N2H"'" maps of 57 low mass cloud cores made at the FCRAO 
14-m antenna equipped with the QUARRY receiver ( [Erickson et al. 1992[ ). The list of 
objects is the same as in BCM98 and it includes all cores already mapped in NH3(1,1). 
The angular resolution of the present observations (54") is 1.5 times finer than the previous 
NH3 study, allowing a more detailed analysis of the morphology and internal motions of 
dense cores. Technical details of the observations made at FCRAO are reported in Sect. 2. 
Results and discussion of this mapping program, including integrated intensity maps, 
column density, size, and mass estimates, velocity gradients, intensity profiles, and the 
variation of line widths across the cores, are shown in Sect. 3. The main conclusions of this 
work are summarized in Sect. 5. 



2. FCRAO Observations 

The J = l—>-0 N2H^ observations were made in March and June 1995 and in March 
1996 at the FCRAO 14 m telescope at New Salem, Massachusetts. We used the 15 element 
QUARRY receiver ( pj^rickson et al. 1992[ ) and the autocorrelation spectrometer with a 



bandwidth of 20 MHz over 1024 channels, giving a spectral resolution of ~20 kHz or 0.063 
km s~^. The beam efficiency (t^b) was 0.51 at 93 GHz ( [Pratap et al. 199^ ) and the main 
beam width at half power (HPBW) was 54". The typical system temperature at 93 GHz 
was ~500 K. The data were acquired in frequency-switching mode with a throw of 8 MHz 
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and calibrated using the ambient load vane method. The rms pointing error, estimated by 
observing Venus and Mars, was ~ 5". 



3. Results and discussion 

3.1. Integrated intensity maps 

We have completed 13 Nyquist sampled maps (25" spacing) and 34 beam sampled maps 
(50" spacing) to intensity levels below half maximum. An additional 10 beam-sampled maps 
are not complete below the half maximum contour. Table [l| reports in columns 2 and 3 the 
coordinates of the (0,0) map position. These do not always correspond to the coordinates 
quoted in BM89 because of recentering of the map to include all the emitting region. The 
peak position is reported in columns 4 and 5 as offsets from the (0,0) coordinates. In 
column 6, the one-sigma level of the noise in the off-line channels (in antenna temperature 
units) is listed. The integrated intensity of the emission at the map peak is in column 7; 
the error on the integrated intensity is o"i = AT^ x ^/iVch x Afj-es, where N^h is the number 
of channels in the integrated intensity (listed in column 8), and Av^es (= 0.063 km s~^) is 
the velocity resolution. Column 9 of Tab. |I| indicates the map grid spacing: beam sampled 
maps (50" spacing) are indicated with a "B", whereas Nyquist sampling (25" spacing) is 
marked by an "N". The size of the mapped area is given in column 10. The association 
with an IRAS source is indicated in column 11. 

A sample of the N2H+(l-0) data is shown in Fig. |l], where the averaged spectrum, 
i.e. the sum of all the spectra inside the half maximum map contour, is presented for 
three starless cores (L1498, L1544, TMC-2) and three cores with stars (L1489, L1228, 
L1251E). The line profiles dramatically change from quiescent starless cores, such as L1498, 
to L1251E, the core with the most complex velocity structure (see Sect. ^.6]). Integrated 



intensity maps are shown in Figure ^. All the maps have the same angular scale and the 
contours are between 20% and 95% of the peak, in steps of 0.15 x the peak intensity. The 
mapped area is shown in the figure (the map type is reported in Tab. |l]). IRAS sources 
associated with cores are marked with stars. The criteria of association are described in 
Jijina, Myers fc Adams (1999) 



3.2. Column density 

The N2H"'"(l-0) line presents hyperfine structure (e.g. Paselli, Myers fc Thaddeus] 



T995t BCM98) and the hfs fitting program in CLASS ( [bbrveille, Guilloteau fc Lucas 1989| ), 
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with the hyperfine frequencies adopted from Caselh et al. 1995, has been used to determine 
LSR velocities (Vlsr); intrinsic hne widths (At>), total optical depths (ttot), and excitation 
temperatures (Tgx). These parameters are listed in Table |^ for the peak spectrum and 
the spectrum averaged inside the half maximum contour. Column 6 of Tab. ^ reports the 
N2H"'" column density (iVtot), which has been calculated following the procedure described 
in [Caselli et al. 2001b[i . In the case of optically thin emission, Tgx = 5 K has been assumed 
to compute A^tot- It is interesting to note that the peak column density averaged over the 
whole sample is (7±5)xl0^^ cm~^, about three times smaller than the averaged column 
density in BCM98 (when the same method is applied). This discrepancy is probably due to 
the smaller beam of BCM98 observations (factor of 2) and suggests the presence of density 
structure in the studied cores (see Sect. [375| ). 

BCM98 found good correlations between N2H"^ and NH3 velocities and line widths, 
but no correlation between column densities of the two species. The authors claimed that 
this result may be due to the different spatial resolutions in the N2H"'" and NH3 studies, 
and to the fact that N2H"'" emission may have a different peak position. In this paper, the 
peak of the N2H"'" emission has been determined and we can then compare the peak column 
densities reported in Tab. |l| with those of NH3. The result is shown in Fig. ^a, where 
only cores with N/aN > 2 have been included. Indeed, the entire sample does not show 
a significant correlation between A^(N2H"'") and A^(NH3). The scatter may be due to the 
significant errors associated with N2H"'" column densities, caused by the large uncertainties 
on the total optical depth probably related to the presence of excitation anomalies (e.g. 
Caselli et al. 1995). However, it is interesting to note that for starless cores the correlation 
between N2H"'" and NH3 column densities is significant (best-fit line in Fig. |]a): 

iVstarless(N2H+) X lO^^ = (2 ± 2) + (0.6 ± 0.3)iVstarless(NH3) X 10^^ Cm-^, (1) 

where the quantities following ± are 1-sigma uncertainties and where the linear correlation 
coefficient is cc = 0.52. The presence of a young stellar object probably affects the chemistry 
in a way that differentiation between the two species starts to be evident. Fig. ^ shows 
that in "NH3-rich starred" cores, the A^(N2H"'")/A^(NH3) column density ratio is smaller 
than in the rest of the sample. This may indicate that the N2H"'"(l-0) line is more optically 
thick and probably more affected by the lower density foreground material, so that a correct 
measurement of r is difficult to obtain without a correct radiative transfer calculation. 

The N2H"'" integrated intensity is plotted versus the "equivalent" NII3 integrated 



^The approximated method used by BCM98, which assumes Jv{Tc^) = Tqx, with Jv{T) the equivalent 
Rayleigh-Jeans temperature, underestimates the N2H+ column density by a factor of about four compared 
to the present calculation. 
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intensity (Ta x Av) in Fig. ^d. Although the associated uncertainties are significantly 
smaller than in the case of column density (see Tab. |^), the dispersion is still large. As in 
the previous case, starless cores show a more significant correlation: 

/stariess(N2H+) = (0.4 ± 0.3) + (2.0 ± 0.7)(rA X At;),tarie..,NH3 Kkms-\ CC = 0.67. (2) 

In Fig. 0c, N2H"'" and NH3 excitation temperatures are plotted. Once again, the correlation 
is strong for starless cores: 

-^ex, starless 

(N2H+) = (1.4±0.8) + (0.4±0.1)Tex ,starless (NH3), (3) 

with cc = 0.79 (dotted line in Fig. ^c), whereas no significant correlation is present in 
cores with stars. Typically, Tex(NH3) > Tex(N2H"'"), which suggests that the critical density 
rZcr(NH3) < ncr(N2H+). In fact, ncr(N2H+) = 2x10^ cm'^ (|Ungerechts et al. 1997|) , about 
one order of mag nitude larger than that of NH3 (2 x 10"^ cm'^; |Swade 19891) . 

Correlations between velocities and linewidths of the N2H~'' and NH3 peak spectra are 
identical to those found in BCM98 and we will not further discuss them. This paper will 
concentrate on the structure and internal motions of individual cores. 



3.3. Sizes and masses 

Sizes of the mapped cores are listed in Tab. ^. The 2D gaussian fitting routine in 
GRAPHIC ( [Buisson et al. 20011 ) has been used to find the position angle, the major and 



minor axis (see columns 2, 3, and 4, respectively). The reported source sizes have been 
corrected for beam size, i.e. we subtracted the gaussian beam size in quadrature for each 
dimension. Aspect ratios (column 5) range between 1.1 and 6.4 with mean ± standard 
deviation 2.0±0.9. The size r in column 6 and 7 is the half power radius, given by the 
geometric mean of the semimajor and semiminor axis. Note that r = R/2, where R is the 
size listed in Table ^ of BM89. The cloud distance is reported in column 8. Comparing the 
size of N2H"'" and NH3 cores with associated stars we find a good correlation (cc = 0.8): 

r,tar(NH3) = (-0.02 ±0.03) + (2.2 ±0.4)rstar(N2H+) PC (4) 

which indicates that the emission of the two molecules have similar morphology, despite the 
different beam sizes. The factor of 2 difference is probably due to differences in the critical 
density ficr of the two tracers since map sizes are already corrected for beam smoothing and 
the resolution ratio is 1.5, less than the typical ratio of radii, 1.8. Indeed, it is interesting 
to note that the above relation pertains to cores with stars, given that for starless cores we 
find: 



?^starless(NH3) = 



(-0.01 ± 0.01) + (0.9 ± 0.3)rstariess(N2H+) pc (cc = 0.7) 



(5) 
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These relations suggest that N2H"'"(l-0) is tracing higher density material than NH3(1,1), 
and that cores with stars are denser and more centrally condensed than starless cores. On 
the other hand, in starless cores, the two lines originate from the same regions. Figure ^ 
shows the distribution of N2H"'" and NH3 core radii in (i) the entire sample, (ii) cores with 
stars, and (iii) starless cores; the two tracers span different size ranges only for cores with 
stars. On average, starless cores have smaller sizes than cores with stars. For 35 cores with 
stars, the N2H+ map radius r has mean ± standard error of the mean (s.e.m.) 0.069±0.005 
pc, while for 19 starless cores, < r > 

starless 

= 0.054±0.005 pc. 

The virial mass of an equivalent uniform density sphere: 

Mvir(M0) = 210 X r(pc) X Au^(km2s-2), (6) 

is listed in column 3 of Tab. ^. The corresponding number density n^i-c is in column 2. Afm 
in eqn. ^ is the FWHM of the molecule of mean mass (2.33 amu, assuming gas with 90% of 
H2 and 10% He): 

^vl = Av' + 8ln2^ X f-l- - , (7) 

where Av is the intrinsic linewidth of the N2H"'" peak spectrum (see Tab. |]), and 'm^2n+ 
(= 29 amu) is the mass of the N2H"'' molecule. Column 4 reports the N2H"'" fractional 
abundance: X(N2H+) = Ar(N2H+)/A^(H2), where A^(H2) = (4/3) x {n^ijl.ll) x r, and 
the factor 1.11 takes into account the difference between nvir, the "virial" number density 
of the molecule of mean mass, and n^[j.{B.2). The N2H"'" column density used for estimating 
X(N2H+) is the peak or the average column density, depending on the associated errors. 
The average value is < X(N2H+) > = (3±l)xl0^i°, close to that found in BCM98. 

Tab. ^ also gives in column 5 the volume density Uex, calculated from the {n^^/n'^^) 
ratio (see eqn. || below), where n'^^ is the critical density (corrected for trapping) of the 
N2H"'"(l-0) line; and in column 6 the "excitation" mass of a uniform core (Mex = (4/3) 
nmnexf^)- The quantity nex/'^cr has been obtained by using the expression for two-level 
statistical equilibrium ( |Genzel 1992| ): 



^ex ^ex -^cb 

n' hE 

k 



where v is the transition frequency, h and k are the constants of Planck and Boltzmann, 
respectively, 

, 1 — exvi—r) , , 

= X ^ (9) 

r 
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and 

f = (^^)(l-exp^-^jj . (10) 

In eqn. ^, r is the optical depth of a "typical" hyperfine component (~ ttot/9). Tex, Tch, 
and Tkin in eqn. || are the excitation, the cosmic background, and the kinetic temperatures, 
respectively. We assumed T^b = 2.7 K, Tkin = 10 K (see BM89), whereas Tex is the 
excitation temperature of the averaged spectrum (see Tab. §). The error on n^x is obtained 
by propagating the error associated with Tex and r into eqn. || (see Appendix 0). 

In Fig. ^ Mex is plotted as a function of M^ir for the whole sample (with the exception 
of cores with (i) M^^/aM^^ or M^y-^jaM^^ < 2, (ii) assumed excitation temperature, and 
(iii) deconvolved size less than the beam size, as in the case of B335). In average, starless 
cores are less massive than cores with stars: < M^ir >star = 9±3 Mq, < M^i^ >stariess = 
3.3±0.4 Mq. However, there is not significant difference in the Mgx/Mvir ratio between the 
two classes of cores (in average < Mex/Mvir >star = 1.4±0.3, and < Mex/Mvir >stariess = 
1.3±0.3). We note that the assumption of a uniform density sphere to estimate the virial 



mass is a very crude one. In Sect. |3.5| we will show that spheroidal cores are consistent with 
density profiles of the form n(r) oc r~°, with < a > ~ 2, although most of the starless cores 
present "central flattening" , at impact parameters b < 5000 AU. If cores were approximated 
with singular isothermal spheres, our virial mass estimates should be reduced by a factor of 
1.6. 



3.4. Velocity Gradients 



Following poodman et al. (1993)| (hereafter GBF93), a least squares fitting of a 



velocity gradient has been performed in all the cores with at least 9 positions with a good 
determination of the LSR velocity Vlsr- In Table ^, the magnitude of the velocity gradient 
Q and its direction {9g, the direction of increasing velocity, measured east of north) are 
reported in column 3 and 4, respectively; the number of velocity points used in the fit is in 
column 2; the product between Q and the core size r, or the typical velocity shift across the 
map, is in column 5; the ratio f3 of rotational kinetic energy to gravitational energy (see 
eqn. (6) in GBF93) is shown in column 6. 

For thirteen of the cores in Tab. |^, the same quantities have been calculated by GBF93 
using NH3 maps (see their Table 2). By comparing the magnitude of the velocity gradient 
calculated from N2H+(l-0) and NH3(1,1) data (^n2H+ and ^nhs, respectively) we found 
no correlation and the majority of the cores (with the exception of L1495 and L1251E) 
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have ^N2H+ > ^NHa, although the scatter is large (the average of the ^n2H+/^nh3 ratio ± 
standard deviation is 1.6±1.0). On the other hand, the correlation between 05(N2H"'") and 
Og(NH3) is significant (cc = 0.7). The slightly larger magnitude of the velocity gradient 
in most of N2H+ cores probably refiects the finer spatial resolution of these observations 
compared to that for NH3. In fact, a larger beam will tend to smooth out Vlsr variations 
across the map. We note, however, that this also implies broader NH3 lines, which are not 



observed (see end of Sect. |3.2| and BCM98). Finally, (3 is similar in N2H"'" and NH3 cores: in 
those cores which are in common in the two samples, the average value of /3 is ~ 0.02 from 
N2H+ data, and ~ 0.03 from NH3 data. 

Tab. 1^ also lists the average value and its standard deviation (columns 7 and 8, 
respectively) of the fit residuals Vlsr(0 " ^t(0 across the core, where Vfnii) is the LSR 
velocity at position (a(z), l3{i)) determined by the least square fit of a linear velocity 
gradient. These quantities are useful to describe the more complex motions in the core and 
will be discussed in section 3.6. 

For 12 cores it has also been possible to determine the magnitude and the direction of 
"local" velocity gradients, i.e. variations of Vlsr in portions of a cloud core. The selected 
cores have at least 12 observed positions where the determination of Vlsr is possible via hfs 
fitting (see Sect. 3.2). In each of these cores, the least square fitting routine to determine Q 
and Qg has been successively applied to adjacent grids of 3x3 points spaced by < 35" or 
< 71" (depending on the map type; see Tab. |l], column 8). Incomplete grids of at least 7 
points have also been included in the computation of "local" velocity gradients. The results 
are shown in Fig. for the 12 selected cores, together with the "global" velocity gradient. 
The arrows point toward increasing Vlsr- From the figure it is evident that many cores 
present internal variations of the magnitude and direction of the velocity gradient (see, 
in particular, L1228 and L1251E), which illustrates the presence of complex motions not 
interpretable as simple solid body rotation of the whole core. 

The velocity gradient maps of Fig. ^ show a wide range of structure which departs from 
the uniform rotation model used by GBF93 to analyse the velocities of the corresponding 
NH3 maps which generally had fewer data points than do the present N2H+ maps. Of the 
9 maps in Fig. |^ having at least 5 velocity gradient vectors, only 2 - L1512 and L1221 - 
have velocity gradient maps which are nearly uniform in both magnitude and direction, 
indicating simple uniform rotation. The rest show patterns which are fairly uniform in 
magnitude but not direction (L183, L1544), uniform in direction but not magnitude (L483), 
or which have significant variation in both magnitude and direction (L1228, L1251E, L43, 
and L1498). The most complex pattern is in L1251E, which shows distinct reversal of 
gradient direction between the East and West parts of the core. The departures from 
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uniform rotation are most pronounced in the cores having associated YSOs and outflows 
(L483, L43, L1228, and L1251E). But substantial departures are also evident in the starless 
cores L183, L1544, and L1498 - which show evidence of contracting motion in CS(2-1) 
( irafaiia et al. 1998| ; [Lee, Myers fc Tafalla 199^1 , 2000), ahhough HCO+(3-2) lines show no 



infall asymmetry toward L1498 ( pregersen fc Evans 20001) . These results are consistent 



with a picture where simple uniform core rotation is rare, where more turbulent motions 
are relatively common, and where such turbulence is associated with cores with YSOs and 
outflows or with starless cores having evidence of inward flows. An interpretation of dense 
core velocity gradients in terms of turbulent motions was recently presented by [Burkert 
Bodenheimer (2000 



We note that cores with complex velocity gradient maps should be poorly fit by a 
simple model of uniform rotation, and inspection of the "normalized" standard deviation 
of the fit residuals ■5<i/LSR-Vflt>/(^ x '^)^ as well as s<i/LSR-Vflt> (see Tab. H), bears out this 
expectation. The most uniform pattern of gradient vectors, for L1512 (Fig. |^), corresponds 
to one of the smallest "normalized" fit residuals, 0.3, while the most complex pattern, for 
L1251E, has one of the largest s ^y^^^_Vfit> I x f)^ 2.4 (and the largest ■s<i/LSH~Vflt>) 0-36 
km s~^). The "normalized" standard deviation relative to L1498 lies between L1512 and 
L1544. This suggests that L1498 is probably in an evolutionary stage later than the "static" 
core L1512 and earlier than the collapsing L1544. In fact, L1498 has been described as an 
extremely quiescent core ( [Lemme et al. 19^ , [Wolkovitch et al. 1997] ) with evidence of slow 



contraction, outer envelope growth and strong chemical differentiations (|Kuiper, Langer 



|fc Velusamy 1996| ; p?afalla et al. 200"l| ). Large values of ■s<VLSR-Vfit>/^ x r (~ 1) are also 



present in TMC-2, an infall candidate with CS asymmetry ( [Lee, Myers fc Tafalla 2000|) . 



3.5. Integrated intensity profiles 

The "standard model" of isolated star formation ([Shu, Adams, fc Lizano 1987[) states 
that cores lose magnetic support by ambipolar diffusion until they become so concentrated 
that the central regions are gravitationally unstable and start to collapse. The collapse of 
the central region deprives the above layers of pressure support and causes them to also 
fall towards the center. In this way, gravitational collapse propagates from the inside out 
and continues until the core runs out of mass or a powerful wind from the central star-disk 
system reverses the collapse and disperses the core. 

Some critical parameters of the "standard model" are unconstrained. The two most 
crucial unknowns are the radial dependence of density and turbulent velocity (alternatively 
of magnetic field), which together determine how gravitational collapse occurs. For example. 
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a change in the density power law from r ^ (isothermal sphere, Shu 1977| ) to r ^ (logotropic 



sphere, [McLaughlin fc Pudritz 1996| ), changes completely the core star forming properties 



( [McLaughlin fc Pudritz 1997] ), and our observations are not fine enough at this point to rule 
out any of these options. In fact, our ignorance of these basic core properties constitutes 
the most serious limitation in our understanding of how stars form in isolated dense cores. 



From submillimeter continuum dust emission, [Ward-Thompson et al. 1994[ found that 



the radial density profiles of pre-stellar cores are significantly different from the singular 
isothermal sphere, and qualitatively consistent with models of magnetically-supported cores 
undergoing ambipolar diffusion (e.g. [Ciolek fc Mouschovias 1995[ ). Typically, the radial 
density profile inferred assuming a constant dust temperature is as flat as p{r) oc r~°, with 
a ~ 0.4-1.2, depending on the core shape, at radii less than ~ 4000 AU, and approaches 



p(r) oc only between ~ 4000 AU and ~ 15000 AU [I 


\ndre, Ward-'i'hompson & Mottej 


|1996 


i Ward-Thompson, Motte, & Andre 1999|; 


Alves et al. 2001 


). However, recent model 


calculations of the dust temperature in pre-stellar cores 


(Zucconi et al. 2001 


; Evans et al. 



2001[) predict a temperature gradient, with a drop from ~ 14 K at the edges to ~ 7 K at the 



centers, which implies more peaked density distributions than in the isothermal case. On 
the other hand, millimeter continuum observations of circumstellar envelopes of low-mass 
protostars are in good agreement with the standard protostellar model of Shu et al. (1987), 
with power-law density gradients such as p(r) oc or r~^'^ ( [Motte fc Andre 2000| ). 

N2H"''(l-0) maps can be used to investigate the column density structure of dense 
cores and make comparisons with results from submillimeter maps. We already noted 
(see Sect. 3.2) that the uncertainty associated with N2H"'" column density is quite large, 
especially for low sensitivity spectra (such as those away from the map peak), because 
of the difficulty in determining an accurate value of the total optical depth of the J=l-0 
transition. Therefore, instead of using N2H"'" column density profiles, we made plots of 
the N2H"*'(l-0) intensity integrated below the seven hyperfine components as a function of 
impact parameter. The use of integrated intensity may be dangerous in those cases where 
the optical depth is large (when the column density is no longer simply proportional to the 
integrated intensity), but we will see that our conclusions are not affected by this problem. 

The integrated intensity profiles of spheroidal cores (with aspect ratio < 2, see Tab. |^) 
have been fitted using two models. Model 1 consists of a spherically symmetric cloud model 
with a density profile p(r) oc r~". The resultant intensity profile vs. impact parameter b (/ 
oc b~^, with p = a — 1), has been convolved with a 2D Gaussian, with FWHM equal to 
the FCRAO beam (54"). For each core, we change the value of p from 0.5 to 2.0, in steps 
of 0.1, and find the best convolved profile and the corresponding a value. Fig. ^ (thin 
curves) and Tab. || (columns 2 and 3) show the results of this procedure. 
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Model 2 considers a spherically symmetric cloud with a radial density profile inferred 
from submillimeter continuum observations of spheroidal cores, with p(r) oc r~^'^ at 
r < Tbreak, and p(r) oc at r > rbreak (see e.g. [Andre, Ward-'l'hompson fc Barsony 2000| ). 
The two profiles have been convolved with the FCRAO beam and joined at rbreak- We run 
several models with different values of impact parameter at the break (febreak from 10" to 
100", in steps of 5"). From the minimization, we determined fcbreak for each core (see 
column 4 of Tab. P). Model 2 profiles are shown in Fig. |^ by the dashed curves. 

The comparison between the values of Model 1 and 2 (columns 3 and 5 of Tab. ^ 
allows one to find the appropriate model density profile for each object. It is interesting 
to note that most (6 out of 9) of the cores with stars are best fitted with single power 
laws. With the exception of Per 6, L1495, and L43, all cores with stars have febreak = 10" 
(the minimum value in Model 2), which is equivalent to having a single slope. Starless 
cores show a different behaviour: 6 out of 8 objects have intensity profiles consistent with 
central flattening at impact parameters less than ~ 5000 AU. This is in good agreement 
with results from submillimeter continuum observations. These results are interpreted 
by |Ward-'l'hompson, Motte, fc Andre 199S as indicating that the cores are probably 



magnetically-supported and evolving through ambipolar diffusion to star formation (e.g. 
Lizano fc Shu 19891 [Basu fc Mouschovias 1995| ). However, caution must be used with 



the above models because, as pointed out by [Andre, Ward-'l'hompson fc Motte 199(j| and 



Andre, Ward-'l'hompson fc Barsony 2000|, they require fairly strong magnetic fields on 



parsec scales (~ 100 /im), difficult to reconcile with available Zeeman measurements (e.g. 
Crutcher 1999[ ). Higher spatial resolution observations are needed to make quantitative 



conclusions on the column and volume density structure of star forming cores and better 
constrain theory. 

In cores with stars, a is typically greater than 2 (Tab. ^). This is probably due to 
an excitation temperature increase at the center caused by the presence of a sufficiently 
luminous protostar. A similar result (greater map "spikiness" of cores with stars over that 
of starless cores) was found by ^ izuno et al. 199^ using H^^CO"'" data. From Fig. ^ we 



also note that the model integrated intensity tends to overestimate data points at large 
impact parameters. This reflects the fast drop in Tex caused by the density drop and it is 
well reproduced by Monte Carlo simulations. It is not related to the "sharp edges" found 
in isolated prestellar cores at radii > 15000 AU ( Pacmann et al. 20001 ). 

To check the effects of optical depth in the observed shallow proflles, in those cores 
with ttot > 10 (L1498, L1489, and L483) we made integrated intensity proflles by using 



•^ttot is the sum of the peak optical depths of the seven hyperfine components 
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the area below the thinnest hyperfine component of N2H+(l-0) having the least optical 
depth (Fi,F = 1,0 1,1; |Uaselh, Myers fc 'I'haddeus 1995| ). In Fig. 0, the dotted curves 
indicate these profiles, normalized so that the peak integrated intesity of the Fi,F = 1,0 
— > 1, 1 component is the same as the peak integrated intensity of all the components. The 
dotted curves in the L1498, L1489, and L483 plots (Fig. 0) closely follow the data points, 
suggesting that optical depth effects are not affecting our conclusions. 

In summary, the radial density profiles of starless cores in Fig. |^ present a consistent 
picture of "central flattening" where the shape of the profile is shallower at small radii than 
at large radii. The N2H"'" integrated intensity profile is modelled by a spherically symmetric 
density law n ~ r~" where a = 1.2 for r < rbreak and a = 2.0 for r > rbreak ~ 0.03 pc. 
Cores with stars are better modelled with single power laws n oc r~" with a > 2. These 
results are the first to show central flattening in a significant sample of molecular line maps. 
In contrast to dust continuum maps, these maps offer the prospect of relating core density 
structure to core turbulence structure, in more detailed studies to be made in the future. 
The agreement between dust and N2H"'" profiles suggests that the molecular gas traced 
by the N2H"'" line is not significantly depleted in relation to the dust, unlike CO and CS 
dUaselh et al. 1999|; [Caselli et al. 2001a|,b; |Bergin et al. 2001|; [lafalla et al. 2001|). 



3.6. Variation of Av across the cores 



Goodman et al. 1998 described a physical picture of star-forming dense cores and 
their environs where the cores are "velocity coherent" regions of nearly constant line width. 
From NH3 maps, Barranco fc Goodman 19TJ^ found that within the interiors of dense cores, 
the line widths are roughly constant and appear to increase at the map edges. Although 
many theories of low mass star formation begin with an isothermal sphere having no 
turbulence (e.g. |!Shu 1977] ), the line width inside the coherent cores is not purely thermal. 
A clearly measurable turbulent component remains even in these "coherent" regions ( Puller 
|fc Myers 1992| ; [Caselli fc Myers 1995|) . The "transition to coherence" may occur because 
of a decrease in the magnetic field's ability to control gas motions in regions of very low 
ionization ( |Mouschovias 1991| ; |Myers 1997| ; [Goodman et al. 1998| ; [Myers 1998|) . 



We have made a similar study with our N2H"'" maps in those cores having a sufficient 
number of high sensitivity spectra across the map to allow us to see systematic variations 
in line width, if they are present. In Fig. ||, integrated intensity contour maps of 9 cores 
selected to have at least 9 spectra with Av/cav > 3 and I/ai > 5 are superposed on grey 
scale maps of line width (light grey indicates narrow lines). From the figure, it is evident 
that starless cores show a spotty pattern of low line width at central positions inside the 
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integrated intensity half-maximum contour, whereas more internal structure is present in 
cores with stars. In fact, L43 ( |Bence et al. 1998|) , L483 ( [I'afalla et al. 2000| , [b'uller fc 



Wootten 2U0D| , [Park et al. TJUUj ) and L1228 ( [1 alalia fc Myers 1997| ) are associated with well 



known protostellar outflows, which may contribute to the broadening of N2H"^ line widths. 
However, in starless cores, the broadest lines are often located at the edges of the map. The 
positions indicated by dots on the maps in Fig. || are those where a good {Av/a/^^ > 3) 
estimate of the intrinsic line width via hfs fitting was possible. 

Another way of looking at core velocity structure is to consider plots similar to those 
shown by [Goodman et al. 1998|, where the intrinsic line width is reported as a function of 



the antenna temperature, which in turn is used as a measure of the distance from the peak of 
the map. Instead of using antenna temperature we consider the impact parameter b, already 
introduced in Sect. 3.5. In practice, at each map point i we associate an "effective radius" 
b{i) by counting the number of positions A^mp with integrated intensity / equal or larger 
than so that b = ^{a x Amp)/?!", where a is the area of the map pixel. Thus we can 
consider in more quantitative fashion how the line width varies with effective map radius. 
We point out that this way of looking at core coherence is equivalent to that described in 
Goodman et al. 1998[ (we tried both methods for three cores and found essentially the same 



results), but this approach shows the dependence on size scale explicitly. 

We first examine how core line widths vary with effective radius by fitting a simple 
linear relation between Av and b. Table ^ (columns 2, 3, and 4) lists the intercept p, the 
slope q, and the correlation coefficient cc of the linear least squares fit to the Av - b data 
for each core (only those cores with at least 9 positions having Av/aAv > 3 and > 5 
are included in the table). From Tab. |^, three classes of cores are recognized: a) cores with 
a "significant" (cc > 0.4) positive Av - b correlation (L1498, L1495, L1524, L1400K, L260); 
b) cores with no significant correlation; and c) cores with a significant negative correlation 
(PER4, B5, TMC-1C2, L1174). Starless cores and cores with stars are found in all three 
classes. Fig. ^ shows some example of the three classes. Evidently the slight increase of 
single-tracer line width with effective map radius discussed by [Goodman et al. 1998| is 
significant in a few cores, but these are not representative of the 22 cores in Tab. 0. 

We next consider how the variation in line width relates to the mean line width, for 
all the usable spectra in the map. Columns 5 and 6 of Tab. |^, report for each core the 
average linewidth < Av > = (1/M) x J^iLi Av{i), with M = number of positions in the 
map where Av has been calculated), and the corresponding standard deviation of the 



sample population s<At,> (= (1/VM- 1) X ^Eiiii^Vi- < Av >)2 ). Note that < Av > 
is different from Av of the average spectrum given in Tab. |[ The corresponding errors are 
reported in Appendix |^. 
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These data can be compared with a simple statistical model of "cells" or "zones" along 
the line of sight to estimate the "coherence length" or length over which the motions are 
correlated along the line of sight ( [Kleiner fc Dickman 19871 ). We assume that each cell 



moves as a coherent unit, with a velocity along the line of sight that follows a Gaussian 
probability distribution with the dispersion a equal to the nonthermal component of the 
overall velocity distribution that we see in the line profile (cx^ = ctnt = Af^/(8Z?T,2) - 
kT/rriohs, where k is the Boltzmann constant, T is the kinetic temperature, and mobs is the 
mass of the observed molecule). Then from basic statistics we can write relations for the 
rms of the line width and the rms of the mean velocity: 

< Ai;nt > .^^s 



a<„> = ' , (12) 

where is the typical number of cells along the line of sight. Figure 10 shows cr<A»;NT> 
versus < Awnt >• In the figure, thin lines represent eqn.( |llD for different values. A 
linear least square fit to the data, taking into account the errors on crA,;NT) gives N = 10, 
with a linear correlation coefficient cc = 0.9. 

Relation ( |T^ ) is investigated in Fig. |TT], where the dispersion of the average velocity 
gradient fit residuals < Vlsr — Vat > is plotted versus the rms of the nonthermal line 



width y< Af^rp >. Vat is the velocity predicted by fitting a first-order gradient, so that 
^LSR — Vfit can be used to analyse higher order structure in the velocity field. In this case, 
a linear least square fit to the data gives = 13, with cc = 0.5. Only cores with at least 
9 data points available to estimate means and standard deviations have been included in 
Fig. |iy and 0. L1251E has been excluded from the two figures because of the complex 
velocity structure clearly seen in Fig. H, suggestive of two adjacent cores rotating in almost 
opposite directions. 



In each of Fig. |10| and |TT], there is a clear tendency for the "dispersion" quantity 
on the ?/-axis to correlate with the "mean" quantity on the x-axis. From this simple 
"cell" model we conclude that N ~ 10. Given that the length along the line of sight to 
which our observations are sensitive is comparable to the map diameter, typically 0.1 pc, 
then the "coherence length" deduced from our data is about 0.01 pc. This size scale is 
comparable to the "cutoff" wavelength below which Alfven waves cannot propagate because 
the neutral-ion collision frequency in the neutral medium becomes comparable to or less 
than the wave frequency ( [McKee fc Zweibel 1995| ; [Goodman et al. 1998| ): 
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where B is the magnetic field strength, Xj is the ionization fraction {rii/n), and n is the 
volume density of the molecule of mean mass (2.33 amu). The choice of i? ~ 10 fxG is based 
on the OH Zeeman measurements in L1544 (a core in our sample) of prutcher fc 'I'roland 



2UUDI At the typical densities of the observed cores (~ 10^ cm ^, see Tab. BI), the ionization 



fraction is about 5x10"*^, if the standard relation between Xi and n{B.2) is used {xi = 1.3 
X 10~^ n{B.2)~^'^, |McKee 198^1) . Therefore, the "transition to coherence" may arise from 
insufficient wave coupling on size scales of ~ 0.01 pc, as (i) proposed by pVIouschovias 1991 



(ii) elaborated by Myers (1997, 1998) as a reason why NH3 line widths are nearly thermal 
and narrower than their surrounding ^'^CO gas (as demonstrated by [Fuller fc Myers 1992| 
and [Caselli fc Myers 1995| ), and (iii) suggested by poodman et al. 1998| as a reason why 



core line widths are "coherent" - nearly constant in the interior of NH3 maps. In general, 
the positive correlations shown in Fig. |TD| and |IT] strongly suggest the presence of internal 
motions, not ascribed to simple solid-body rotation and probably due to turbulence, which 
broaden lines and contribute to line width and LSR velocity dispersion across the cores. 

To further quantify the "coherence" of core line widths, we ask whether the "local" 
dispersion of line widths at radius near b, sav, increases with b in N2H"'" cores. To answer 
this question, the standard deviation of the average Av inside all the contour levels between 
20% and 100% of the integrated intensity peak has been calculated (in steps determined 
by the following requirements: i) consecutive levels are associated with different impact 
parameters; ii) each level contains at least 5 positions). As with previous analysis, we only 
considered those positions where Av/a^v > 3 and I/cxi > 5. We then performed linear 
least square fits to the s^v^b data in each core and the results are listed in columns 7, 8, 
and 9 of Tab. ^. Note that all the cores (with the exception of L1495, L1536, and L1512) 
show a strong positive correlation between the two quantities, indicating that dispersion 
is indeed increasing with projected radius. We tested whether the general increase in 
with increasing b could be due to the corresponding decrease in signal-to-noise (S/N) ratio 
with increasing b. We raised the threshold (//o"i)niin from 5 to 10 and repeated the linear 
lest-square fits as in Tab. with smaller, lower-noise data sets. For the 12 surviving cores 
with enough data points (at least 9), the tendency for sad to increase with b is evident in 5 
starless cores - L1544, L183, TMC-ICS, TMC-1NH3, TMC-2 - and 5 cores with stars - 
L1489, L43, L1228, L1221, L1251E. The remaining two cores show a negative (L1512) and 
a null (L483) sa?) — b correlation. 



We note that Caselli et al. 200rb| have recently found a different Xi-n(H2) relation for the L1544 core: Xi 



— 5.2xl0~^ n(H2)~'~''^^ , which, at n(H2) = 10^ cm~'^, implies an electron fraction of 8xl0~^, six times lower 
than that deduced from the "standard" relation. With this new value of Xi, the "typical" cutoff wavelength 
is about 0.001 pc. 
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This analysis of 22 low-mass N2H"'" cores indicates that the brightest part of a 
core map has nearly constant line width, typically 0.3-0.4 km s~^, while the surrounding 
positions with fainter emission have line widths with much greater variation. This confirms 
the conclusions of Barranco fc Goodman llJ^ based on 4 NH3 core maps. However, our 
data show no significant increase in N2H"'" line width with effective map radius for the 
typical core, in contrast to the result of Parranco fc Goodman 1998. 



3.7. Map structure 

The N2H"'" maps presented in Fig. |^ are generally similar in shape and orientation 
to the NH3 maps of BM89 and |Ladd, Myers fc Goodman 1994| . This similarity is easily 
understood, since N2H"'" and NH3 are closely related chemically (N2 is the precursor 
molecule for both species). In cores with stars, the factor of ~ 2 smaller size of the N2H"'" 
maps compared to its NH3 counterpart probably reflects both the higher critical density 
(~ 2x10^ cm"'^ for N2H+ compared to ~ 3x10^ cm~'^ for NH3) and the finer angular 
resolution (~ 54" instead of ~ 83") and sampling (25" or 50" instead of 60") used in the 
N2H+ observations. 

To quantify the complexity of the projected map structure, we have counted the 
number of map "peaks" and "elongations" in each of the 59 N2H+ maps in this paper and 
in each of the 47 NH3 maps in BM89 and in Ladd, Myers fc Goodman 1994| . Here a "peak" 



is a local maximum of integrated intensity which exceeds its surrounding valley by at least 
3-sigma, and an "elongation" is an extension from the position of a peak by at least one 
beam, which does not terminate in a new peak. We find that the N2H"*" and NH3 maps 
have similar proportions of simple and complex structure. The fraction of " simple" maps 
with single peaks and no elongation is 0.64 for NH3 and 0.70 for N2H+ . The fraction of 
maps with at least two peaks is 0.12 for NH3 and 0.07 for N2H"'" . Because of small number 
statistics, these NH3 and N2H"'" fractions are not significantly different. One might expect 
that more peaks per map would be seen in the finer-resolution N2H+ observations than in 
the coarser-resolution NH3 observations, but this is not the case. We interpret this result 
as arising from the fact that the N2H+ maps have both finer resolution and smaller spatial 
extent, each by about the same factor of 1.5. Thus some N2H"'" maps show two peaks where 
the NH3 map shows one, but other N2H+ maps do not extend far enough to sense the 
second peak seen in the NH3 map. Although relatively few N2H"'" maps have more than one 
local maximum (14 of 57), we note that most cases with double peaks have peak-to-peak 
separation of only 1-2 FWHM beam diameters. Therefore these core maps are nearly as 
"clumpy" as the map resolution allows. 
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3.8. Overview of core parameters 

An overview of tlie parameters determined in previous sections is presented in Table ^ 
for cores with stars and starless cores separately. Although the dispersion in large, cores 
associated with young stellar objects are typically more turbulent (non-thermal line widths 
are about 1.5 times larger), have larger sizes (factor of ~ 1.4), and are more massive (factor 
of ~ 2), than starless cores. Other quantities (Tex, A^tot(N2H+ ), aspect ratio, Q, jS, 
have very similar values in the two classes of cores, suggesting that there is not a definite 
separation between them, at least from the analysis of (1^0) lines at the present 

angular resolution. 



4. Conclusions 

We have mapped 57 low mass cores in the rotational transition J = l—>-0 of the 
molecular ion N2H"'", using the FCRAO antenna. This extensive mapping survey has 
allowed us to study physical properties of dense cores with an angular resolution about 1.5 
times finer than previous ammonia maps from BM89. The main conclusions of this work 
are summarized below. 

1. The excitation temperature of the N2H"'" (1-0) line is typically ~ 5 K, indicating that 
N2H+ lines are subthermally excited. The peak N2H+ column density averaged over the 
whole sample is A^(N2H"'" ) ~ 7x10^^ cm~^, about two orders of magnitude less than A^(NH3 
) (BM89). There is a positive correlation between N2H+ and NH3 column densities and 
excitation temperatures in starless cores, whereas in cores with stars the scatter is large and 
no significant correlations are found. Although this may partially be due to the difficulty 
in estimating the N2H+(l-0) total optical depth, the lack of correlation in cores with stars 
suggests a different chemical evolution of NH3 and N2H"'". However, the good correlations 
between LSR velocities and line widths in the entire sample, indicates that the two tracers 
generally originate from the same regions. 

2. The mean aspect ratios of the mapped sources is ~ 2. Starless cores have about the 
same linear sizes than those found with NH3 (1,1) maps. On the other hand, cores with 
stars have r(NH3 ) ~ 2 x r(N2H+ ). This gives evidence that cores associated with young 
stellar objects are more centrally concentrated than starless cores and that N2H+ (1-0) 
traces denser gas than NH3 (1,1). 

3. In average, starless cores have virial and "excitation" masses M^^j. ~ M^x ~ 3 Mq, and 
are less massive than cores with stars {M^[^ ~ M^x ~ 8 Mq). The M^x/M^ij. ratio averaged 
over the whole sample is ~ 1. 
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4. Typical values of velocity gradient magnitudes are ~ 2 km/s/pc, both in cores with stars 
and starless cores. If the gradient represents rotation, the ratio (3 of rotational energy to 
gravitational energy ranges between ~ 10~^ and 0.07 so that rotation is not significant in 
the support of the core. Maps of "local" velocity gradients reveal the presence of complex 
internal motions that deviate strongly from a simple model of solid body rotation of the 
whole core. 

5. Six out of nine spheroidal starless cores present central flattening in the integrated 
intensity profile. This is consistent with a spherically symmetric density law n{r) ~ r~° 
where a = 1.2 for r < rbreak and a = 2 for r > rbreak) with rbreak ~ 0.03 pc. Cores with 
stars are better modelled with single power law density profiles with a ^ 2. These results 
are in qualitative agreement with submillimeter continuum observations, suggesting that 
N2H"'" is not significantly depleted inside dense cores (unlike CO and CS). 

6. Most N2H"'" cores are "coherent" in having more uniform line widths in their bright 
interior than in their faint periphery, as seen earlier in NH3 observations. The fluctuations 
in line width also increase significantly with mean line width. However this sample shows 
no significant tendency for the line widths themselves to increase with map radius - a few 
cores have positive and negative trends, while most have no significant trend. For 20 of 26 
cores, the standard deviation of the average line width, s^y, increases with b, indicating 
that core line widths vary more with increasing radius as previously found by [Barranco 



&; Goodman 1998| and poodman et al. 1998| using ammonia data. Yet, line widths Av 



positively correlate with the impact parameter b in only 5 sources (L1498, L1495, L1524, 
L1400K, and L260). Four sources present negative correlations (PER4, B5, TMC-1C2, 
L1174). The remaining 17 cores do not show a significant Av — b correlation. 

7. The "coherence length" deduced from our data is about 0.01 pc, comparable to the 
cutoff wavelength below which Alfven waves cannot propagate. Thus, the "transition to 
coherence" may arise from a decay of turbulence in the innermost parts of the cores, due to 
insufficient wave coupling, on size scales of ~ 0.01 pc. 

8. Although N2H+ maps have finer angular resolution than NH3 maps, they do not show 
a more complex structure. The majority (70%) of the cores in our sample have "simple" 
N2H"*" maps, with single peaks and no elongation. Most cases with double-peaks have 
peak-to-peak separation of only 1-2 FWHM beam diameters. 
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A. Error expressions 

The error on the volume density rzex (see Sect. p73|) is found by propagating the error 
in eqn.(P): 



k fkin(fkin - Tcb) T4exp(-/i«/(fcTox)) 1 - e- 



T2 



3t 



ncx3r e-3^(3r + l)-l 
1 - T 



1/2 



(Al) 



The errors on the average hne width < Av > and the corresponding sample standard 
deviation s^v have been calculated with the following expressions (from the propagation of 



error : 



C" <Ad> 



S<Ai;> 



N 



S<Av>iN-l) 



(A2) 
(A3) 
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Table 1. N2H+(l-0) map results: peak position {Aa, A6) and integrated intensity (/) 
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RA(1950)'^ 


Dec(1950)^ 


Aa 
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A,5 
W 
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(K) 
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B 
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n 
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B 
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B 


21.4 


y 


PER7 


03 29 39.5 


30 49 50 


0.00 


0.00 
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B 
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0.04 
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B 
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0.00 
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0.19 
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95 


B 
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04 00 38.0 


56 47 59 
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B 
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0.16 
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98 


N 
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25 02 13 


0.42 
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0.16 


1.20±0.10 


88 


N 
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Table 1 — Continued 
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0.00 
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0.88 


0.80 


0.18 


1.91±0.12 
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0.17 
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N 


42.8 


y 
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23 23 32.2 


74 01 45 


0.86 


-0.04 


0.15 


1.64±0.10 
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B 
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y 



^Coordinates of the (0,0) map position; they do not always correspond to the coordinates quoted in BM89. 

"^Integrated intensity at the peak position; tj/ = AT,^ x v^/Vch x A?;rcs, where AT^ ^.j^^^ [K] is the 1 a level of the 
noise in the off-line channels, iVch is the number of channels in the integrated intensity, and Avr^s [km s~^] is the velocity 
resolution (=0.063 km s~^). 

'^Type of map: B = beam sampling; N = Nyquist sampling. 

■^Mapped area. 

°The half maximum contour extends over the mapped area. 

^A second component is present in the West direction; it has not been included in the integrated intensity estimate. 

^There is another peak in the mapped area which belongs to another core (whose half maximum contour extends outside 
the mapped area). 

^IRAS association, from Jijina et al. 1999. 
'Number of channels in integrated area. 
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Table 2. N2H+(l-0) map results: multicomponent (hfs) fit to the peak and averaged*^ 

spectrum. 



Core 




Av 




rp C 

ex 


Ntotx 10-12 




(km s~^) 


(km s~^) 




(K) 


(cm-2) 


PER4-Ad 


7.60±0.02 


0.58±0.06 


0.1 


5.0 


4.7±0.6 




7.60±0.02 


0.47±0.05 


4±2 


5±1 


5±3 


PER4-B 


7.61±0.01 


0.36±0.03 


0.1 


5.0 


4.5±0.5 




7.565±0.007 


0.36±0.01 


0.1 


5.0 


3.5±0.1 


PER4-C 


7.38±0.02 


0.33±0.04 


0.1 


5.0 


3.0±0.4 




7.41±0.02 


0.38±0.04 


0.1 


5.0 


2.5±0.4 


PER5 


8.221±0.008 


0.33±0.02 


0.1 


5.0 


5.5±0.5 




8.203±0.006 


0.29±0.02 


5±2 


5±1 


5±2 


PER6 


5.87±0.01 


0.44±0.03 


0.1 


5.0 


4.8±0.4 




5.880±0.008 


0.38±0.02 


6±2 


4.1±0.5 


6±2 


PERT 


6.81±0.02 


0.33±0.04 


15±7 


3.8±0.6 


11±5 




6.82±0.02 


0.45±0.05 


5±3 


3.9±0.8 


5±3 


PER9 


6.91±0.02 


0.41±0.04 


0.1 


5.0 


3.4±0.4 




6.82±0.01 


0.42±0.03 


0.1 


5.0 


2.9±0.2 


B5 


10.31±0.02 


0.47±0.03 


0.1 


5.0 


6.0±0.5 




10.25±0.01 


0.39±0.02 


4±2 


4.7±0.9 


5±2 


L1389 


-4.65±0.03 


0.49±0.05 


0.1 


5.0 


3.2±0.4 




-4.60±0.01 


0.33±0.04 


7±3 


3.9±0.7 


5±3 


L1489 


6.80±0.01 


0.28±0.02 


17±7 


4.3±0.8 


13±5 




6.783±0.004 


0.277±0.009 


7±1 


4.7±0.4 


6±1 


L1498 


7.833±0.007 


0.18±0.02 


18±9 


4.1±0.8 


8ib4 




7.840±0.002 


0.191±0.005 


llibl 


3.9±0.2 


4.7±0.6 


L1495 


6.832±0.007 


0.22±0.02 


4±3 


7±3 


5±3 




6.824±0.004 


0.27±0.01 


4±1 


6±1 


4±1 


LMOOG^ 












B217 


7.02±0.01 


0.34±0.03 


4±3 


6±3 


6±4 




7.001±0.005 


0.33±0.01 


4±1 


5.1±0.8 


5±2 


L1524 


6.36±0.01 


0.26±0.03 


15±8 


4.1±0.9 


10±5 




6.339±0.008 


0.35±0.02 


9±2 


3.7±0.3 


7±2 


L1400K 


3.28±0.01 


0.19±0.02 


11±6 


3.6±0.6 


4±2 




3.300±0.006 


0.24±0.01 


6±2 


3.6±0.3 


3±1 


TMC-2A 


5.917±0.006 


0.22±0.02 


16±4 


4.7±0.7 


11±3 




5.935±0.006 


0.27±0.01 


9±2 


4.2±0.4 


7±2 


TMC-2 


6.26±0.02 


0.40±0.03 


0.1 


5.0 


4.9±0.5 




6.197±0.005 


0.36±0.01 


5±1 


4.3±0.4 


5±1 
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Table 2 — Continued 
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IS.V 


TTOT*^ 


7^ c 

ex 


Ntotx 10-12 




(km s"-*-) 


(km s~^) 




(K) 


(cm-2) 


L1536 


5.53ib0.01 


0.30±0.02 


r\ 1 

0.1 


5.0 


3.9ib0.4 




5.649±0.003 


r\ r^/^o 1 r\ r\f\c\ 

0.263±0.009 


6±1 


4.2ib0.3 


3.9±0.8 


L1534 


A f\ 1 r\ f\'\ 


0.36zb0.02 


r\ 1 

0.1 


5.0 


4.6±0.4 




6.304ib0.006 


c\ or\ 1 (\ c\-\ 


1 1 

3ibl 


r 1 1 

5ibl 


A 1 C\ 

4±2 


L1527 


f f\c\ 1 c\ c\-\ 

5.90±0.01 


0.29±0.03 


13±7 


3.9±0.7 


8ib5 




5.922±0.008 


0.34±0.02 


9±2 


3.6ib0.3 


6±2 


TMC-INH3 


v f^fr* 1 r\ r\r\r\ 

5.956ib0.009 


0.36ib0.02 


llib3 


4.6±0.5 


llib3 




v fw v 1 r\ r\c\ A 

5.955±0.004 


0.39±0.01 


5.6±0.8 


4.6±0.3 


6.2±0.9 




5.27±0.02 


0.29±0.03 


14±D 


3.9±0.7 


9±4 




5.219ib0.006 


c\ c\c\ 1 (\ r\i 

0.29ib0.01 


8±2 


3.8ib0.3 


r 1 1 

5ibl 


TMC-IC 


v rtT 1 c\ c\c\ 

5.27±0.02 


0.27±0.04 


18±10 


3.8ib0.8 


11±6 




5.257±0.004 


0.25±0.01 


12±2 


3.6±0.2 


6±1 


TMC-ICS 


5.89±0.01 


0.38±0.03 


8±3 


4.3±0.8 


8±3 




5.853±0.005 


C\ A 'I 1 /"»/"» 1 


A i 1 

4±1 


4.3±0.4 


5=bl 


L1517B 


5.80±0.01 


0.27±0.02 


r\ 1 

0.1 


5.0 


3.0±0.3 




5.830ib0.008 


0.22±t).t)2 


r\ 1 A 

9±4 


4.0ib0.6 


5±2 


L1512 


7.108±0.007 


0.18±0.02 


5±3 


5±2 


3±2 




7.088±0.002 


0.195ib0.006 


7±1 


4.2ib0.3 


3.6±0.6 


LI 544 


7.169ib0.008 


r\ 1 1 1 
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10ib2 


4.5±0.5 
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7.162±0.003 


0.307±0.006 


8.1±0.9 


4.0±0.2 


6.0±0.7 


T 1 r" 00 A I 

LI 582 A 


1 C\ C\C\ 1 C\ C\C\ 

10.20±0.02 


0.43±0.05 


r\ 1 

0.1 


5.0 


1.9ib0.3 


B35 


11.d9±0.04 


0.61±0.08 


r\ 1 A 
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3.6±0.5 


12ib6 




11.8d±0.04 
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r\ 1 
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2.74±0.01 
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0.34±0.02 
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L183 


2.44±0.01 


0.25±0.02 


22±7 


4.3±0.7 


14±5 
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0.303±0.009 


8±1 


4.1±0.2 


6.1±0.9 


L1681B-A 


3.62±0.02 


0.40±0.03 


0.1 


5.0 


5.5±0.6 




3.69±0.02 


0.45±0.04 


0.1 


5.0 


3.8±0.4 


L1681B-B 


4.12±0.01 


0.36±0.03 


0.1 


5.0 


5.7±0.7 




4.13±0.01 


0.39±0.04 


0.1 


5.0 


3.9±0.5 


L1696A 


3.37±0.01 


0.25±0.03 


0.1 


5.0 


3.9±0.6 




3.40±0.01 


0.30±0.02 


0.1 


5.0 


2.9±0.3 


L43 


0.678±0.009 


0.36±0.02 


6±2 


7±1 


11±4 
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Table 2 — Continued 



Core 




Av 


TTOT^ 


c 

^ ex 


Ntotx 10-12 




(km s~^) 


(km s~^) 




(K) 


(cm-2) 




0.692ib0.003 


0.396ib0.007 


5.9±0.5 


5.4ib0.3 


8.6±0.8 


L260 


3.503ib0.008 


0.20±0.02 


0.1 


5.0 


2.2ib0.3 




3.480ib0.005 


0.19zb0.01 


llzb3 


3.4ib0.2 


A 1 1 

4±1 


T 1 r" of 

L158 


3.91ib0.01 


0.23ib0.04 


0.1 


5.0 


2.7ib0.5 


L234A 


2.95±0.01 


0.21ib0.02 


0.1 


5.0 


2.3±0.3 




2.936±0.007 


0.23±0.02 


0.1 


5.0 


1.9±0.2 


L234Eg 














3.04±0.02 


0.30±0.04 


12±7 


3.2±0.3 


7±4 


T O 

L63 


5.78ib0.01 


0.21±0.03 


13ib6 


r 1 1 

5±1 


8±4 




5.780±0.007 


0.27±0.02 


9±2 


4.3ib0.5 


D±2 


B68 


3.35±0.02 


0.27±0.04 


0.1 


5.0 


1.9ib0.4 




3.38ib0.02 


0.33±0.05 


0.1 


5.0 


1.4ib0.3 


L483 


5.53±0.01 


0.59±0.03 


16±3 


4.6±0.4 


27=b5 




5.432±0.0l)4 


0.495±0.008 


11.6±0.7 


4.5±0.1 


1d±1 


B133 


12.1±t).2 


0.9±0.3 


C\ 1 

0.1 


5.0 


C\ 1 1 

2±1 


L778 


9.98 ±0.02 


0.44±0.04 


6ib3 


A A \ C\ C\ 

4.4±0.9 


7±4 




9.944±0.009 


0.40±0.02 


5±2 


4.4ib0.7 


5ib2 


B335 


8.36±0.01 


0.39±0.03 


6±3 


4.5±0.9 


6±3 




8.345ib0.008 


0.40±0.02 


7ib2 


A C\ \ C\ A 

4.0±0.4 


D±2 


L1152 


2.70ib0.02 


0.4d±0.04 


r\ 1 

0.1 


5.0 


4.2±0.D 




2.60±0.02 


0.48±0.03 


0.1 


5.0 


2.8±0.2 


L1155C 


2.69ib0.02 


0.36ib0.04 


7ib4 


3.6±0.5 


5±3 




2.70±0.01 


0.33ib0.03 


6±3 


3.5±0.4 


4ib2 


L1082C 


-2.53±0.02 


0.42±0.04 


0.1 


5.0 


3.8±0.5 




-2.55±0.01 


0.33±0.03 


16±6 


3.4±0.3 


10±4 


L1082A 


-2.13±0.02 


0.35±0.04 


0.1 


5.0 


1.9±0.3 




-2.19±0.02 


0.46±0.04 


0.1 


5.0 


1.7±0.2 


L1228 


-8.06±0.02 


0.61±0.04 


9±2 


4.9±0.6 


18±4 




-8.041±0.008 


0.68±0.02 


5.9±0.8 


4.5±0.3 


11±2 


L1174 


2.67±0.06 


1.2±0.2 


0.1 


5.0 


7±1 




2.66±0.02 


1.05±0.06 


0.1 


5.0 


4.8±0.4 


BERN48 


-7.35±0.03 


0.54±0.06 


0.1 


5.0 


5.4±0.8 




-7.34±0.01 


0.46±0.04 


4±2 


5±1 


5±3 



-30- 



Table 2 — Continued 



Core 


^LSR 




b 

TTOT 




NtotXlO 




(km s ) 


(km s ) 




(K) 


(cm ^) 


L1172A 


2.91±0.04 


54±0 09 


8±5 


3.8±0.7 


10±6 




2 88±n 02 


50±0 04 


8±3 


3.7±0.4 


9±3 


B361 


2 78±n 05 


70±0 09 

\J • 1 \J—l—\Jm\JtJ 


0.1 


5.0 


3.5±0.6 




2 fi4±0 04 


0.9±0.2 


5±3 


3.3±0.3 


9±5 


LIOSIC^ 












L1031B 


4 18±0 05 


1.6±0.2 


0.1 


5.0 


7.0±0.9 




o.yi±U.Uo 


1 1 -i-n 1 


n 1 
U.i 


o.U 


4.0±U.O 


L1221 


-4.42±0.02 


0.67±0.05 


6±2 


6±1 


15±5 




-4.430±0.009 


0.69±0.02 


6.9±0.9 


4.5±0.3 


13±2 


L1251A 


-3.93±0.02 


0.45±0.05 


0.1 


5.0 


4.1±0.6 




-3.93±0.01 


0.36±0.04 


8±4 


3.8±0.6 


6±3 


L1251C 


-4.71±0.01 


0.29±0.03 


12±5 


4.3±0.9 


9zb4 




-4.75±0.02 


0.55±0.06 


4±3 


3.9±0.9 


6±4 


L1251E 


-3.93±0.03 


1.35±0.08 


0.1 


5.0 


11.8±0.8 




-3.872±0.009 


1.50±0.02 


O.liO.l 


5.0 


7.8±0.2 


L1262 


4.11±0.01 


0.39±0.05 


18±10 


3.8±0.8 


15±9 




4.06±0.01 


0.43±0.02 


8±2 


3.8±0.3 


8±2 



^The averaged spectrum is obtained by adding together all the spectra inside the half 
maximum contour. 

^ttot is the sum of the peak optical depth of the seven hyperfine components. 
'^Excitation temperature calculated by assuming a main beam efhciency 7/6=0.51. If ttot < 
1, Tex = 5 K has been assumed (see text). 

'^The first row refers to the peak spectrum, whereas the second row refers to the averaged 
spectrum. 

^No detection. 

^Compact source: only one spectrum inside the half maximum contour. 
^Individual spectra cannot be hfs-fitted because of low S/N. 

^Only the peak spectrum can be hfs-fitted because of low S/N in the other map spectra. 
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Table 3. Angular and linear size of N2H"'"(l-0) cores 



Core^ 


u A b 
FA 


Major*^ 


Minor 


Aspect 




r 


U 




(deg) 


(arcmin) 


(arcmin) 


Ratio 


(arcmin) 




(pc) 


PER4-A' 


45 


3.4 


0.5 


6.4 


0.7 


0.069 


350 


PER4-Bf 


98 


3.9 


0.8 


4.6 


0.9 


0.092 


350 


PER4-C 


22 


2.0 


1.3 


1.6 


0.8 


0.08 


350 


PER5 


59 


2.2 


1.3 


1.7 


0.8 


0.09 


350 


PER6 


64 


3.1 


1.5 


2.0 


1.1 


0.11 


350 


PER7 


40 


1.7 


1.0 


1.6 


0.7 


0.07 


350 


PER9 


71 


1.8 


1.1 


1.7 


0.7 


0.07 


350 


B5 


''98 


3.9 


0.9 


4.3 


0.9 


0.095 


350 


L1389^ 


22 


1.2 


0.5 


2.5 


0.4 


0.065 


600 


L1489 


135 


1.6 


1.1 


1.4 


0.7 


0.027 


140 


L1498 


121 


2.6 


2.0 


1.4 


1.1 


0.046 


140 


L1495 


115 


2.9 


1.9 


1.5 


1.2 


0.048 


140 


B217 


58 


2.7 


1.6 


1.7 


1.0 


0.042 


140 


L1524 


71 


5.3 


1.8 


2.9 


1.6 


0.064 


140 


L1400K 


35 


2.9 


1.8 


1.6 


1.2 


0.057 


170 


TMC-2A 


89 


2.0 


1.4 


1.4 


0.8 


0.034 


140 


TMC-2 


53 


5.3 


3.8 


1.4 


2.2 


0.091 


140 


L1536 


162 


5.5 


3.2 


1.7 


2.1 


0.085 


140 


L1534 


147 


5.7 


2.0 


2.8 


1.7 


0.069 


140 


L1527 


19 


2.5 


1.9 


1.3 


1.1 


0.044 


140 


L1517B 


165 


1.6 


1.4 


1.1 


0.7 


0.030 


140 


L1512 


114 


2.3 


1.6 


1.5 


1.0 


0.039 


140 


L1544 


125 


2.0 


1.2 


1.7 


0.8 


0.032 


140 


L1582Af 


19 


1.2 


0.5 


2.5 


0.4 


0.046 


400 


B35 


f3 


1.8 


0.9 


2.1 


0.6 


0.072 


400 


LI 34 A 


121 


3.3 


2.0 


1.7 


1.3 


0.060 


160 


L1681B-A 


146 


2.3 


1.0 


2.2 


0.8 


0.036 


160 


L1681B-B 


162 


2.4 


1.1 


2.2 


0.8 


0.039 


160 










J. .o 


n Q 




1 fin 


L43 


172 


3.0 


1.6 


1.9 


1.1 


0.050 


160 


L260 


38 


3.7 


2.1 


1.8 


1.4 


0.065 


160 


L158 


^112 


0.9 


0.8 


1.2 


0.4 


0.020 


160 


L234A 


45 


3.6 


1.2 


3.0 


1.0 


0.048 


160 


L63 


94 


3.7 


1.9 


1.9 


1.3 


0.062 


160 


B68 


97 


2.1 


1.1 


1.9 


0.8 


0.044 


200 


L483 


47 


1.5 


1.2 


1.2 


0.7 


0.040 


200 


B133 


30 


1.3 


1.0 


1.3 


0.6 


0.101 


600 


L778 


148 


2.0 


1.6 


1.2 


0.9 


0.109 


420 


B335f 


78 


0.7 


0.3 


2.1 


0.2 


0.018 


250 


L1152 


*46 


3.4 


0.9 


3.6 


0.9 


0.115 


440 
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Table 3 — Continued 



Core'* 




Major"^ 


Minor 


Aspect 




r 






(deg) 


(arcmin) 


(arcmin) 


Ratio 


(arcmin) 


(pc) 


(pc) 


L1155C 


80 


1.7 


1.1 


1.6 


0.7 


0.088 


440 


L1082C 


^45 


1.2 


0.8 


1.6 


0.5 


0.062 


440 


L1082A 


81 


1.6 


1.2 


1.3 


0.7 


0.092 


440 


L1228 


95 


1.5 


1.0 


1.5 


0.6 


0.055 


300 


L1174 


165 


1.3 


1.0 


1.3 


0.6 


0.072 


440 


BERN48f 


8 


1.3 


0.5 


2.9 


0.4 


0.023 


200 


L1172A 


28 


1.7 


1.5 


1.2 


0.8 


0.100 


440 


B361 


156 


2.7 


1.6 


1.7 


1.1 


0.108 


350 


L1031B 


133 


1.5 


1.0 


1.5 


0.6 


0.164 


900 


L1221 


126 


1.8 


1.0 


1.8 


0.7 


0.039 


200 


L1251A 


58 


2.7 


1.4 


2.0 


1.0 


0.056 


200 


L1251C 


45 


1.8 


1.3 


1.4 


0.7 


0.043 


200 


L1251E 


1 


5.2 


1.8 


2.9 


1.5 


0.088 


200 


L1262 


155 


2.1 


1.4 


1.5 


0.9 


0.050 


200 



''Six cores do not appear in this table because their half maximum contours extend beyond the mapped 
area: TMC-INH3, TMC-1C2, TMC-IC, TMC-ICS, L183, and L234E. In the cases of L1534, L134A, 
and L260 it was possible to fit the core with a 2D Gaussian because only a small fraction of the half 
maximum contour lies beyond the mapped area. 

'^The position angle PA is defined as the angle in a clockwise direction from the positive right ascension 
axis. 

'^The major and minor axes have been corrected for beam size. 

'^r is the half-power radius, 0.5 times the geometric mean of the major and minor axis. Note that r = 
R/2, where R is the size listed in Tab. | of BM89. 

"Distance references are listed in BM89, Ladd et al. (1994), GBF93, Jijina et al. (1999). 

^These cores have deconvolved sizes similar to the beam size, so their small sizes are less certain than 
the sizes of the larger sources because of beam subtraction. 
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Table 4. Volume Density and Mass 



Core 


(10^ cm-3) 




X(N2H+)'' 
(10-10) 


(10^ cm-3) 




PER4-A 


10±1 


8±1 


1.9±0.4 


15±11 


15±11 


PER4-B 


3.2±0.3 


6.0±0.5 


4.1±0.6 


23 


49 


PER4-C 


3.9±0.4 


4.9±0.5 


2.6±0.5 


23 


32 


PER5 


3.2±0.2 


5.3±0.2 


5.0±0.5 


18±8 


30±14 


PER6 


2.8±0.2 


8.6±0.6 


4.3±0.5 


8±3 


28±10 


PER7 


5.6±0.5 


4.1±0.4 


8±4 


4±2 


4±1 


PER9 


6.0±0.6 


5.3±0.5 


2.2±0.3 


23 


23 


B5 


3.9±0.3 


8.1±0.6 


4.4±0.5 


15±8 


35±19 


L1389 


9±1 


5.7±0.7 


1.5±0.3 


7±4 


6±4 


L1489 


32±1 


1.47±0.06 


4±2 


13±3 


0.7±0.1 


L1498 


8.8±0.3 


2.08±0.07 


5±3 


5.7±0.7 


1.4±0.2 


L1495 


8.8±0.3 


2.33±0.09 


3±2 


28±14 


8±4 


B217 


15±1 


2.6±0.2 


3±2 


20±8 


4±2 


L1524 


5.4±0.3 


3.3±0.2 


7±4 


5±1 


3.1±0.8 


L1400K 


5.9±0.2 


2.61±0.09 


3±2 


5±2 


2.3±0.7 


TMC-2A 


17.7±0.7 


1.64±0.06 


5±1 


8±2 


0.8±0.2 


TMC-2 


3.6±0.3 


6.5±0.5 


4.0±0.5 


11±3 


20±5 


L1536 


3.3±0.1 


4.9±0.2 


3.8±0.4 


10±2 


14±3 


L1534 


5.7±0.3 


4.5±0.2 


3.2±0.3 


23±13 


19±10 


L1527 


12.0±0.8 


2.5±0.2 


4±2 


4±1 


0.9±0.3 


TMC-1NH3 








13±2 




TMC-1C2 








6±1 




TMC-IC 








3.8±0.6 




TMC-ICS 








11±3 




L1517B 


25±1 


1.61±0.07 


l.liO.l 


7±3 


0.5±0.2 


L1512 


12.5±0.4 


1.74±0.06 


2±1 


9±2 


1.3±0.2 


L1544 


23.8±0.5 


1.87±0.04 


3.1±0.7 


7±1 


0.60±0.08 


L1582A 


16±2 


3.5±0.4 


0.7±0.1 


23 


8 


B35 


9±2 


8±1 


5±2 


4±2 


4±2 


L134A 


6.6±0.4 


3.3±0.2 


1.4±0.2 


23 


12 


L183 








8±1 




L1681B-A 


24±2 


2.6±0.2 


1.8±0.2 


23 


3 


L1681B-B 


18±1 


2.5±0.2 


2.2±0.3 


23 


4 


L1696A 


13.3±0.8 


2.1±0.1 


2.0±0.3 


23 


4 


L43 


11.0±0.5 


3.3±0.2 


6±2 


22±3 


6.8±0.9 


L260 


4.6±0.2 


3.0±0.1 


2.0±0.3 


3.0±0.7 


2.0±0.5 


L158 


50±4 


1.00±0.08 


0.7±0.2 


23 


0.6 


L234A 


8.7±0.3 


2.27±0.08 


1.5±0.2 


23 


7 


L234E 








1.9±0.8 




L63 


5.2±0.3 


2.9±0.2 


7±4 


9±3 


5±2 
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Table 4 — Continued 



Core 




Mvir'' 


X(N2H+)i^ 








(10^ cm^3) 


(Mq) 


(10-1") 


(10^ cm-3) 


(Mq) 


B68 


12±1 


2.3±0.2 


1.0±0.2 


23 


5 


L483 


29±2 


4.5±0.3 


6±1 


9.3±0.5 


1.55±0.08 


B133 


8±5 


19±12 


0.7±0.6 


23 


65 


L778 


2.8±0.3 


8.6±0.8 


6±3 


12±5 


39±18 


B335 


89±6 


1.27±0.09 


1.1±0.5 


8±2 


0.21±0.06 


L1152 


2.6±0.2 


9.5±0.9 


3.7±0.6 


23 


93 


L1155C 


3.6±0.3 


5.7±0.5 


4±3 


4±2 


8±4 


L1082C 


8.2±0.8 


4.7±0.4 


2.0±0.3 


2.5±0.7 


1.7±0.5 


L1082A 


3.2±0.3 


5.9±0.5 


1.8±0.3 


23 


47 


L1228 


16±1 


6.4±0.6 


5±1 


12±2 


5.4±0.9 


L1174 


26±7 


23±6 


0.9±0.3 


23 


24 


BERN48 


79±11 


2.3±0.3 


0.8±0.2 


15±10 


0.6±0.4 


L1172A 


4.1±0.8 


10±2 


6±4 


5±2 


13±5 


B361 


5±1 


15±3 


1.7±0.4 


3±2 


10±5 


L1031B 


9±2 


95±18 


1.3±0.3 


23 


276 


L1221 


36±4 


5.2±0.6 


3±1 


12±2 


1.8±0.3 


L1251A 


11±1 


4.5±0.5 


1.8±0.4 


6±3 


3±1 


L1251C 


12.4±0.8 


2.4±0.2 


4±2 


8±4 


1.7±0.8 


L1251E 


23±2 


37±4 


1.6±0.2 


23 


39 


L1262 


12±1 


3.5±0.4 


7±4 


6±1 


1.9±0.5 



'^TT-vir and Mvir is the virial volume density and mass, respectively. Data are 
not reported for those cores where the size cannot be determined. 

''Fractional abundance of N2H+ (X(N2H+) = iV(N2H+)/iV(H2)) calculated 
from nvir and assuming a uniform sphere with A^(H2) = 4/3 nvir/1.1 r (the 
factor 1.1 is to convert n to n(H2)). 

■^Volume density and mass coming from the density to critical density (ricr 
= 2x10^ cm-3; [Ungerechts et al. 1997|, ApJ, 482, 245) ratio, calculated by 



using equation (43) in Genzel 1992 , Values with no associated errors imply 
an assumed Tex value (— 5 K). The excitation temperature and the optical 
depth of the "averaged" or the "peak" spectrum (see Table 2) have been used, 
whichever has the smallest error. In the calculation, Tkin = 10 K and Tbb — 
2.7 K. 
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Table 5. Results of Gradient Fitting 







a 


^6 


^ X V 


P 


^ t'LSR •'fit 


*<l^LSR-^tit> 






( l^TTl / G / T^l^ 1 
1 JVlll f O f ^Kj 1 


CHps- F, of N") 


(km /s^l 






Ckm /si 


PER4 


10 


0.64±0.05 


-17±6 






-O.OliO.Ol 


0.07±0.01 


B5 


11 


0.86±0.08 


84±2 


0.082 


9.2 


0.015±0.007 


0.07±0.01 


L1489 


19 


0.7±0.1 


110±12 


0.018 


0.74 


0.004±0.004 


0.02±0.01 


L1498 


26 


0.5±0.1 


9±10 


0.024 


1.4 


-0.002±0.003 


0.013±0.003 


L1495 


9 


0.9±0.2 


87±8 


0.045 


4.5 


0.004±0.006 


0.046±0.008 


B217 


11 


2.2±0.2 


36±4 


0.093 


16 


0.024±0.008 


0.060±0.008 


L1524 


9 


1.8±0.3 


22±11 


0.11 


29 


0.012±0.009 


0.03±0.01 


L1400K 


10 


1.8±0.1 


62±4 


0.10 


27 


0.000±0.005 


0.019±0.006 


TMC-2 


24 


0.7±0.1 


86±7 


0.067 


6.6 


-0.018±0.009 


0.093±0.009 


LI 536 


21 


2.11±0.07 


-7±1 


0.18 


66 


-0.002±0.004 


0.041±0.004 


L1534 


15 


2.5±0.2 


25±3 


0.17 


53 


-0.011±0.007 


0.083±0.008 


TMC-1NH3 


16 


5.98±0.08 


35.1±0.7 






-0.008±0.005 


0.159±0.006 


TMC-1C2 


14 


1.25±0.09 


94±6 






0.004±0.005 


0.036±0.007 


TMC-IC 


15 


0.8±0.1 


-136±10 






0.004±0.005 


0.040±0.005 


TMC-ICS 


17 


2.7±0.2 


44±3 






-0.003±0.006 


0.030±0.008 


L1512 


28 


1.41±0.07 


166±3 


0.054 


7.7 


-0.002±0.002 


0.017±0.003 


L1544 


32 


l.OiO.l 


-171±5 


0.032 


2.0 


0.011±0.003 


0.033±0.003 


L183 


35 


1.19±0.08 


-55±3 






-0.006±0.003 


0.055±0.005 


L43 


42 


2.72±0.06 


127±2 


0.13 


33 


-0.022±0.003 


0.045±0.004 


L260 


9 


O.litO.l 


86±94 


0.0078 


0.11 


-0.022±0.007 


0.04±0.01 


L483 


44 


2.38±0.06 


63±2 


0.096 


9.5 


0.022±0.004 


0.063±0.004 


B335 


10 


0.1±0.3 


-58±187 


0.0014 


0.005 


-0.005±0.008 


0.021±0.009 


L1228 


43 


2.2±0.1 


-64±2 


0.12 


15 


0.048±0.009 


O.lOitO.Ol 


L1174 


13 


4.0±0.2 


-158±4 


0.29 


30 


0.01±0.02 


0.08d=0.03 


L1221 


22 


3.2±0.3 


-138±5 


0.13 


14 


0.02±0.01 


0.07±0.01 


L1251E 


90 


1.66±0.06 


-129±3 


0.15 


5.8 


-0.175±0.007 


0.355±0.006 



'^(3 is the ratio of rotational kinetic energy to gravitational energy (see equation 6 of Goodman et al. 1993). 
Assuming po = m X rivir, with m = 2.33 amu, and rivir = volume density from Table 4, /3 = 4.86x10^ < 6^ > /n^\i 
with Q in km/s/pc. 

''Average value (and relative standard deviation) of the fit residuals VlsrC*) " ^t(i) across the core, where Vflt(«) 
is the LSR velocity at position (a(«), 5{i)) determined by the least square fit of a velocity gradient. 
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Table 6. Integrated intensity profiles 



Core 


a 


xl 


^brcak 


xl 






(Xl02) 


(pc) 


(Xl02) 


Cores with stars 


PER6 


1.9 


5.8 


0.051 


3.7 


L1489 


2.2 


0.65 


<0.007'* 


8.7 


L1495 


1.9 


3.4 


0.024 


1.3 


L43 


1.8 


31 


0.019 


23 


L483 


2.2 


23 


<0.010 


55 


L1082C 


2.5 


0.24 


<0.021 


9.4 


L1228 


2.8 


23 


<0.015 


220 


L1174 


2.2 


0.61 


<0.021 


8.8 


L1221 


2.1 


22 


<0.010 


28 






Starless 


cores 




L1498 


1.8 


2.4 


0.014 


3.0 


B217 


1.9 


1.9 


0.020 


0.86 


L1400K 


1.9 


0.32 


0.021 


0.33 


TMC-2 


1.7 


15 


0.041 


6.5 


L1536 


1.6 


16 


0.037 


3.0 


L1512 


1.9 


6.4 


0.020 


3.3 


L1544 


1.9 


2.1 


0.014 


0.96 


L63 


1.8 


13 


0.031 


4.8 



^The smallest is obtained with the 

smallest 6break (10")) indicating that the 

profile is consistent with a single power law 
(see text for details). 
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Table 7. Variation of Av across the cores 



Core 






cc 




S<Av> 


pic 




cc' 




(km s-l) 


(km pc"'^) 




(km s-^) 


(km s-^) 


(km s-l) 


(km pc~^) 








n Q-t-n A 




n 'i7-\-c\ no 
u.o 1 mu.uz 


u. lumu.uo 


n 9-t-n 1 


n A-t-O Q 
-u.^mu.y 


n Qi 


-DO 


U-^OIIIU-UO 


1 s-t-n 'i 
-i-omu.o 




U-OOZCU-Ui 




u.u^mu.uo 


n A-t-O A 


n 01 
u . y i 


L1489 


U.24±U.U2 


1.4±0.7 


0.30 


U.298±0.U08 


O.Ud±O.U2 


n no 1 n n 1 

0.U3±0.U1 


0.2it0.4 


0.23 


L1498 


0.19±0.01 


1.2±0.4 


0.41 


0.235±0.006 


0.038±0.009 


0.008±0.010 


0.8±0.3 


0.81 


L1495 


0.21±0.03 


2.0±0.8 


0.53 


0.29±0.01 


0.05±0.01 


0.05±0.06 


Oil 


0.11 


B217 


0.33d=0.03 


-0.5±0.8 


-0.14 


0.33±0.01 


0.07±0.02 


-0.01±0.05 


1±1 


0.68 


L1524 


0.24±0.04 


2±1 


0.52 


0.34±0.02 


0.07±0.02 


O.litO.l 


-1±2 


-0.95 


L1400K 


0.18±0.03 


1.3±0.6 


0.80 


0.24±0.01 


0.02±0.01 


0.04±0.04 


-0.3±0.8 


-0.92 


TMC-2 


0.35±0.03 


-0.7±0.5 


-0.14 


0.40±0.01 


0.11±0.02 


0.01±0.02 


1.4±0.4 


0.88 


L1536 


0.25±0.02 


-0.2±0.3 


-0.08 


0.248±0.007 


0.044±0.008 


0.04±0.01 


0.1±0.2 


0.20 


L1534 


0.38±0.03 


-0.7±0.6 


-0.25 


0.36±0.01 


0.05±0.01 


0.00±0.02 


0.7±0.4 


0.91 


TMC-1NH3 


0.40±0.02 


-1.3±0.4 


-0.30 


0.39±0.01 


0.13±0.02 


-0.08±0.03 


3.2±0.5 


0.79 


TMC-1C2 


0.34±0.03 


-1.5±0.6 


-0.60 


0.282±0.009 


0.04±0.01 


0.01±0.03 


0.6±0.6 


0.54 


TMC-IC 


0.26±0.03 


O.litO.6 


0.03 


0.29±0.01 


0.05±0.01 


0.03±0.03 


0.5±0.5 


0.54 


TMC-ICS 


0.38±0.03 


0.2±0.6 


0.06 


0.40±0.01 


0.07±0.02 


0.03±0.03 


0.5±0.6 


0.62 


L1512 


0.17±0.01 


0.8±0.4 


0.29 


0.202±0.005 


0.033±0.006 


0.019±0.007 


0.3±0.2 


0.70 


L1544 


0.30±0.01 


-0.1±0.4 


-0.03 


0.303±0.006 


0.043±0.008 


0.017±0.007 


0.6±0.2 


0.82 


L183 


0.23it0.01 


0.7±0.3 


0.29 


0.282±0.008 


0.06±0.01 


-0.005±0.009 


0.9±0.2 


0.85 


L43 


0.27±0.01 


1.7±0.3 


0.30 


0.394±0.007 


O.llitO.Ol 


0.031±0.006 


1.2±0.1 


0.84 


L260 


0.18±0.03 


1.3±0.8 


0.51 


0.24±0.01 


0.05±0.02 


-0.01±0.08 


1±2 


0.59 


L483 


0.40±0.01 


-0.4±0.3 


-0.07 


0.45±0.01 


0.15±0.02 


0.067±0.009 


1.0±0.2 


0.65 


B335 


0.43±0.04 


-1±1 


-0.27 


0.39±0.02 


0.07±0.02 


0.05±0.08 


0±2 


0.43 


L1228 


0.65±0.03 


-0.7±0.4 


-0.16 


0.68±0.02 


0.19±0.02 


0.12±0.02 


0.7±0.3 


0.61 


L1174 


1.2±0.1 


-9±2 


-0.73 


0.81±0.04 


0.36±0.07 


0.4±0.3 


0±3 


0.49 


L1221 


0.70±0.04 


-2±1 


-0.24 


0.69±0.02 


0.15d=0.04 


0.00±0.04 


3±1 


0.90 


L1251E 


0.90±0.02 


-3.0±0.3 


-0.31 


0.93±0.02 


0.34±0.02 


0.25±0.01 


l.OibO.l 


0.80 



is the intercept and q is the slope of the best-fit Hnear Aii — b relation in each core (see text). 
''Average linewidth and corresponding standard deviation. 

'^p' is the intercept and q' is the slope of the sav - b relation in each core (see text). 
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Table 8. Statistics on cores with and without stars 







Cores with St 


ars 




Starless Cores 




Mean 


Standard 


Number 


Mean 


Standard 


Number 


Parameter 




Deviation 


of Cores 




Deviation 


of Cores 


AvtiT (km s~^) 


0.5 


0.3 


35 


0.3 


0.1 


25 


Tex (K) 


5 


1 


15 


4.4 


0.8 


14 


^tot(N2H+) (1012 cm-2) 


8 


5 


35 


6 


3 


25 


r (pc) 


0.07 


0.03 


22 


0.05 


0.02 


13 


Aspect ratio 


1.9 


0.7 


35 


2 


1 


19 


rivir (10^ cm-3) 


2 


2 


34 


1.2 


0.7 


19 


M^i, (Mq) 


9 


16 


34 


3 


2 


19 


nex (IQS cm-3) 


0.9 


0.7 


21 


0.8 


0.4 


17 


Mex (Mq) 


8 


11 


21 


3 


6 


17 


X(N2H+) (10-10) 


3 


2 


34 


2 


1 


18 


Gradient (km/s/pc) 


2 


1 


14 


2 


1 


12 


/3 


0.02 


0.02 


13 


0.02 


0.02 


7 



'^Adnx is the non-thermal part of the line width; Ari^^ = Ad^|^_^ - 8ln(2) x (feT/mobs)- 
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Fig. 1. — Average N2H+(l-0) spectra of selected cores, obtained by summing all the spectra 
inside the half maximum map contour. The spectra are in antenna temperature units. 
Three starless cores (L1498, L1544, TMC-2) and three cores with stars, indicated by the 
symbol "*", are displayed. This sub-sample show the variation in the N2H"'" profile from 
quiescent (L1498) to more "turbulent" starless cores (L1544, TMC-2), and from relatively 
quiescent cores with stars (L1489), to cores with young stellar objects driving powerful 
outflows (L1228), to cores with complex internal structure (L1215E). 

Fig. 2. — Maps of the N2H"'"(l-0) intensity integrated over the seven hyperfine components. 
The angular scale is the same for each core. The contours and the grey scale mark the 20%, 
35%, 50%, 65%, 80%, and 95% of the map peak, reported in Table 1 (see column 7). The 
thick contour is the half maximum (50%) level, which defines the core size. Small circles 
are observed positions and the stars indicate the location of the associated infrared source 
detected by IRAS. The FCRAO half power beam width (HPBW) is shown in the map of 
Per 4. 

Fig. 3. — Correlations between N2H+ and NH3 (from BM89) properties at the map peaks 
for the cores with distance D < 200 pc: (a) column density; (b) integrated intensity; (c) 
excitation temperature. Best-fit fines are referred to starless cores (black dots). Cores with 
stars (grey dots) do not show significant correlations between N2H+ and NH3 properties (see 
text). 

Fig. 4. — Distribution of core radii mapped in N2H"'" (shaded histogram) and NH3 (thin 
lines), for the entire sample (top), cores with stars (center), and starless cores {bottom). 
Cores associated with young stellar objects tend to have smaller sizes in N2H+ than in NH3 
maps, suggesting higher central densities than starless cores. 

Fig. 5. — The "excitation" mass M^x, calculated assuming a spherical core with constant 
density n^x (sec text), as a function of the virial mass M^ir for starless cores (empty circles) 
and cores with stars (filled circles). Cores with Mex/cTMex or Mvir/cxM^i, < 2 are not reported 
in the figure. 

Fig. 6. — N2H"''(l-0) integrated intensity maps of those cores where "local" velocity gradients 
have been calculated. The grey-scale levels represent the 30%, 50%, 70%, and 90% of the 
map peak. Small dots mark the position of observed spectra where the determination of Vlsr 
from hfs fit has been possible (sec Sect. 3.2). The white arrows show the magnitude and the 
direction of the velocity gradient calculated by applying the least square fitting routine to 
the grid of positions centered on the corresponding arrow. The black arrow in the bottom 
of each panel represents the total velocity gradient listed in Table 5 (the magnitude of the 
white arrows is in units of the total gradient). 
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Fig. 7. — N2H"'"(l-0) integrated intensity as a function of impact parameter b. Symbols are 
different for starless cores (empty circles) and cores associated with stars (filled circles). Thin 
(dashed) curves are Model 1 (Model 2) best-fit profiles (see text). Dotted profiles have been 
obtained by using the normalized integrated intensity of the "thin" hyperfine component 
(Fi,F = 1,0 1,1). Most of the starless cores present a "shallow" structure at 6 < ~ 0.03 
pc, in agreement with results from dust continuum emission maps. 

Fig. 8. — Maps of line width (grey scale) overlapped with integrated intensity maps (contours; 
levels are 30, 50, 70, and 90% of the peak). Grey contours range from Av^am to < Av^aax in 
steps of 2 < aAv >, where < ctav > is the mean line width error in the selected positions. 
Grey areas enclose all the points with Av values between two adjacent grey contours. The 
dots mark the positions which have been used in the Av maps (i.e. where Av/a/}^^ > 3 
and //cT/ > 5). Values of Avmin, ^^max, and 2 < a/^^ > (in km s~^) in each core are the 
following: i) 0.15, 0.29, 0.05 in L1498; ii) 0.21, 0.53, 0.11 in TMC-2; iii) 0.13, 0.24, 0.05 in 
L1512; iv) 0.19, 0.39, 0.07 in L1544; v) 0.18, 0.38, 0.08 in L183; vi) 0.24, 0.70, 0.08 in L43; 
vii) 0.24, 0.85, 0.11 in L483; viii) 0.32, 1.28, 0.21 in L1228; ix) 0.37, 1.91, 0.27 in L1251E. 

Fig. 9. — Intrinsic N2H"*"(l-0) line width Av as a function of impact parameter b for selected 
cores with positive (L1498, L1495), null (TMC-2, L1228), and negative (TMC-1C2, B5) 
correlations (see text for details). Cores with stars are marked with filled symbols whereas 
starless cores have emtpy symbols. 

Fig. 10. — The variance of the average nonthermal line width of a core as a function of 
< Awnt >• Empty circles represent starless cores, whereas filled circles are cores with stars. 
Lines indicate how dispersion is expected to increase with < Af nt > in a simple model of 
"cells" along the line of sight (see text). The number of cells is indicated. In both classes of 
cores, dispersion is increasing with increasing < Awnt >, following models with N ~ 10. 

Fig. 11. — Dispersion of the average gradient fit residual (Tab. ^ in a core vs. the average 
nonthermal line width (Tab. |^. The "cell" model is indicated by thin lines for number of 
cells equal to 1, 3, 10, and 30, as in Fig. |10|. A linear least square fit to the data gives N 
= 13. These positive correlations strongly suggest the presence of turbulent motions which 
cause line width broadening and contribute to line width and LSR velocity dispersion across 
the cores. 
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